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Summary. — In these lecture notes, the role of strangeness in relativistic as- 
trophysics of compact stars is addressed. The appearance of strange particles, as 
hyperons, kaons, and strange quarks, in the core of compact stars is examined and 
common features as well as differences are presented. Impacts on the global prop- 
erties of compact stars and signals of the presence of exotic matter are outlined for 
the various strange phases which can appear in the interior at high densities. 

PACS 26.60,+c - Nuclear matter aspects of neutron stars. 

PACS 97 . 60 . Jd - Neutron stars. 

PACS 21 . 80 . +a - Hypernuclei. 

PACS 12.38.Mh - Quark-gluon plasma. 
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1. Introduction 

Neutron stars are the final cndpoint in the evolution of stars more massive than eight 
solar masses. They are born in a spectacular explosive event, a core-collapse supernova, 
which can outshine an entire galaxy in its brightness. A surprisingly good estimate about 
the characteristic masses and radii can be made by adopting the argument of Landau [1] . 
The delicate balance between gravitational energy and Fermi energy can only be achieved 
up to a certain critical number of nucleons (fermions) which is N s=s (Mp/ttin) 3 w 10 57 , 
where M P stands for the Planck mass and m,/v for the nucleon mass. The corresponding 
maximum mass for neutron stars amounts to M max w M 3 /m 2 N k, 1.8M q . This mass 
estimate is exactly in the region of maximum masses presently discussed in the literature 
with much more refined theoretical models of neutron star matter. There are excellent 
modern textbooks on compact stars available, which discuss the physics of neutron stars 
and its various phases in great detail. We refer the interested reader to the books by 
Norman Glendenning [2] and by Fridolin Weber [3] for a thorough treatment of this 
exciting field in the interplay of physics and astrophysics. Most recently, a new textbook 
on compact stars appeared which updates and complements the existing ones [4]. A 
recent review article on the relation between the nuclear equation of state and neutron 
stars can be found in [5] , a nice popular article on neutron stars in [6] . 

In these lecture notes, we discuss the global properties of neutron stars, masses and 
radii, and how those change for different compositions in their interior. In particular, we 
address the appearance of new and exotic phases in dense neutron star matter and how 
they can modify the mass-radius diagram for compact stars. In section 2, the baryons 
of the SU(3) octet are treated together in a chiral SU(3)xSU(3) effective Lagrangian 
which is motivated from the approximate symmetries of the underlying theory of strong 
interactions, quantum chromodynamics (QCD). When the properties of hyperons are 
tied to the known experimental data on hypernuclei, hyperons are present in neutron 
star matter above twice normal nuclear matter density affecting the properties of massive 
neutron star configurations. In section 3, strange mesons, the K~ , and its presence in 
neutron star matter by forming a Bosc-Einstcin condensate are addressed. A relativistic 
field-theoretic approach is utilised to outline the basic features of kaon condensation 
for compact stars. Ultimately, hadrons have to be described in terms of quark degrees 
of freedom at high densities. The phase transition from the chirally broken hadronic 
phase to the approximately chirally restored quark matter phase is examined within 
perturbative QCD and the MIT bag model to illustrate the main features of a strong 
first order phase transition in compact star matter in section 4. Characteristics of the 
mass-radius relation are highlighted which signal the presence of exotic matter in the 
core of compact stars. 

2. Hyperons in Neutron Stars 

QCD with massless quarks is chirally symmetric, which means that all left-handed and 
right-handed quarks decouple. This statement is actually true for all vector interactions. 
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Left- and right-handed quarks are defined by the relations 



(1) 



Ql = ^(i -75)9 
1, 

QR = ^(1 +75j9 



(3,0) 



(0,3) 



Splitting a spinor * to left and right-handed components one gets 



(2) 



§7^* = {L + Rh^{L + R) = L^A»L + Rj^R 



The mass term for quarks violates chiral symmetry as 



(3) 



m** = m(L + R)(L + R) = m(LR + RL) 



and quarks with different chirality mix with each other. The explicit breaking of chiral 
symmetry is small, however, as m U: d ~ 10 MeV and also m s ~ 100 MeV is much smaller 
than the nucleon mass to at w 1 GeV so that chiral symmetry is a useful tool. QCD 
has a complex and nontrivial structure, quark and gluon condensates are present in 
the vacuum, in particular the light quark condensate a, the strange quark condensate 
£ and the gluon condensate \. The nonvanishing vacuum expectation values of these 
condensates actually generate most of the masses of the hadrons (except for the pseudo- 
Goldstone bosons). 

In constructing an effective chiral Lagrangian for the description of neutron star mat- 
ter, one has to consider composite quark fields for the meson and baryon fields. First, 
consider the spin zero mesons. Assuming that they are s-wave bound states, then the 
only spinlcss objects we can form are 



The combinations q L qh and q R qn vanish, since the left and right chiral subspaces are 
orthogonal to each other. The resulting representation in chiral SU(3) x SU(3) sym- 
metry is then (3,3*) and (3*, 3), respectively. The antiparticles belong to the conjugate 
representation. Hence, nonets of pscudoscalar and scalar particles have to be considered. 
For the vector mesons, one has to construct vector-like quantities out of the quark fields 
qi and q R . Again assuming s-wave bound states, the only vectors which can be formed 
are 

(5) IlI^l q R iv<iR 



(4) 



U1R 



This suggests assigning the vector and axial vector mesons to the representation (3 x 3*,0) 
(0,3 x 3*) = (8,1) 8 (1,8), an octet and a singlet state. 
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The meson fields can be grouped conveniently into matrices under flavour SU(3), here 
for the scalar fields and pseudoscalar fields 

8 8 

(6) ^(g L A a te + q L X a lb q R ) = ^(£ a A a + ^ a A a ) = S + ill = M 

a=0 a=0 

and correspondingly for the vector and axial vector fields as well as for the baryon fields. 
Out of all these fields one constructs a Lagrangian which obeys the chiral SU(3) symmetry 

(7) C = £kin + £BM + £BV + £vcc + Co + CsB + Aep , 

with the usual kinetic terms for baryons (£kin), spin-0 fields (Co), vector mesons (£ VC c) 
and lcptons, electrons and muons (C\ cp ). The baryons couple with Yukawa-type inter- 
actions to the spin-0 mesons (£bm) and to spin-1 mesons (£bv)- An explicit chiral 
symmetry breaking term is introduced also (£sb)- For homogeneous matter and in the 
mean-field approximation, derivatives of the boson fields vanish and the fields with un- 
natural parity (pseudoscalars and axial vector fields) vanish. One arrives at the following 
expressions for the various terms in the Lagrangian: 

Cbm + Cbv = - [m* + giulou a + 5i07o0° + 9NpJoT 3 po] ipi , 

i 

r 1 2X 2 2,1 2^i2 , 2 2, 4/ 4 , oi4 , c 2 2 4\ 

C vcc = -mi— u' + -m% — r + -— ™^+g|(w 4 + 20 4 + 6wV+P ) , 
z Xo z Xo z Xo 

Co = -\k x 2 (<j 2 + C 2 ) + h(<T 2 + C 2 ) 2 + fc 2 (y + C 4 ) + h X a\ 

,4 1 4l xV 8. a\ 
-kiX ~ tX In —a + o ln ■ 



Csb = — — 



2 

X 



4 A X 4 o 3 a 2 Co ' 
ml.U<J + (V2m 2 K f K - -^=mlj^)C 



(8) A C p = - mi]tpi . 



Nonvanishing vacuum expectation values are generated by spontaneous breaking of chiral 
symmetry, so there appears a light quark condensate and a strange quark condensate. 
In addition, effects from the nonvanishing gluon condensate are taken into account by 
introducing a scalar \ field, which is a chiral singlet. These condensates generate the 
hadron masses. So hadron masses are not additional input parameters in this effective 
field theory and come out automatically by considering chiral SU(3) symmetry, contrary 
to e.g. the standard relativistic mean-field theory (there is for example no explicit mass 
term for baryons in the Lagrangian!). The expressions for the baryon masses read for 
example 



m 



N 



m - \gos{^OS - 1)(V2C - o-) 
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Fig. 1. - The hadron masses as generated by the vacuum expectation values of the chiral SU(3) 
model in comparison to the experimental data. 



m* A = m - -.9o 8 (a os - 1){V2( - a) 
2 

m| = m + -.9o 8 (a os - 1){V2( - a) 
(9) ml = m + -gf) 8 {2aos + 1)(V2C - a) . 

with mo = gQi(y/2cr + £)/V3. The generated hadron mass spectrum is compared to the 
experimental data in Fig. 1. 

For determining the necessary input to the Tolman-Oppcnheimer-Volkoff (TOV) equa- 
tion, the relation between pressure and energy density has to be calculated. The ther- 
modynamic grandcanonical potential 



can be derived by standard methods. The particle energy of the baryons depends now 



6 



J. SCHAFFNER— BlELICH, S. SCHRAMM, and H. STOCKER 



on the expectation values of the meson fields in the medium 
(10) E;(k i ) = yfkZ + mt* 

so that the baryons acquire an effective mass m* and an effective chemical potential 

(f f ) Mi = E *{ k F,i) + SW^o + 9i<t><f>o + 9i P hiPo , 

which is shifted by the vector potentials. All thermodynamic quantities, as the number 
density n and the energy density e, can be extracted from the thermodynamic potential 
via 

,„„. ^ dp \ 

(12) P= ~V ' ni = ~dp~ ' e = -P + 2^W«i • 

The coupling constants of the scalar mesons to the baryons are fixed by determining 
the baryon masses in vacuum. The vector coupling constant are automatically given by 
SU(3) symmetry relation: 

cn\ o ^ v gs<t> V% 

(lo) gAui — g^Lu — £g~.ui — g.9Afw = 2g 08 ; 9A<f> = = = ~^-gNu> ■ 

which are actually the SU(6) symmetry relations known from the quark model (ideal 
mixing is assumed which is a very good approximation for the vector meson nonet). 

The chiral effective model gives a good description of nuclear matter as well as of 
the properties of nuclei. The nuclear matter properties are a binding energy of Eb/A = 
— 16 MeV at n = 0.15 fm~ 3 with an effective mass of m* N /mN = 0.61, a compression 
modulus of K — 276 MeV, and an asymmetry term of a sym = 40.4 MeV. Hyperons are 
automatically included by adopting consistently SU(3) symmetry for the chiral effective 
Lagrangian from the beginning. The computed single-particle energy levels of hypernuclei 
are depicted in Fig. 2 and compared with the experimental hypernuclear data. One sees 
an overall agreement with the data, just reflecting the fact that A hypernuclei are well 
described by a potential depth of the A at saturation density of about —30 MeV. The 
situation for the other hyperons, £ and S, is far less clear. The £ hypernuclear potential 
is likely to be repulsive, while there is experimental evidence that the S hypernuclear 
potential is attractive, although significantly less in comparison to the one for A hyperons. 
In the chiral effective model, the hyperon potential are fixed already by the parameters 
and the SU(3) symmetry relations. The hyperon potential for the S comes out to be 
only barely repulsive in the chiral model, which will be important for the composition of 
neutron star matter. 

Hyperons can appear in neutron star matter by virtue of the conditions of /3-equilibrium 
[9, 10]. The timescale for weak interaction rates is of the order of 10~ 10 s (decay timescale 
for hyperons) to 10~ 8 s (decay timescale for kaons) in vacuum. In matter, similar rates 
are expected, as for example hyperons in hypernuclei have similar to just slightly smaller 
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Fig. 2. - The single-particle energy levels of various hypernuclei for the SU(3) chiral model 
compared to the experimental data (taken from [7]). 



lifetimes compared to the hyperon lifetime in vacuum. Neutron stars are known to 
be quite old, some of them have characteristic ages close to the age of the universe 
of 10 10 years, plenty of time to reach equilibrium for weak interactions, /3-equilibrium. 
Hence, A hyperons can appear in dense neutron star matter if their effective energy in 
the medium reaches the baryochemical potential of neutrons at some baryon density n: 
EK 71 ) — MA — M«- The relations between the chemical potentials of all particles are 
determined by the conserved charges of the particles, which for neutron star matter are 
just baryon number and charge (not strangeness, contrary to e.g. heavy-ion reactions). 
Hence, the chemical potential of all particles can be fixed by 

(14) m = Bi ■ hb + Qi ■ 

where B t and Qi are the baryon number and the charge of the particle i, respectively. 
The effective energy for hyperons increases less rapidly than the one for nucleons, due 
to the SU(3) symmetry of the vector coupling constant of the vector potential, which 
dominates the high-density behaviour. Therefore, modern calculations of the composition 




Fig. 3. - The composition of neutron star matter as a function of density for the chiral SU(3) 
model (upper plot) and the relativistic mean-field model using the parameter set TM1 (lower 
plot). Hyperons appear around twice normal nuclear matter density (taken from [8]). 
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of neutron star matter predict that hyperons appear around 2n , where n stands for 
the normal nuclear matter density. Fig. 3 shows the composition as calculated from 
the chiral effective model presented above. Neutrons are the main component at low 
densities. Protons and electrons appear then with equal amounts, so as to conserve 
charge neutrality. Then at about 2.5no first the £~ then the A hyperons are present in 
the matter, reaching fractions of around 10% at larger densities. The £~ appears before 
the A despite its heavier mass as negatively charged particles are favoured in neutron star 
matter, so as to balance the positive charge of the protons. The asymmetry energy of 
nucleons drives the system to isospin-neutral matter, i.e. equal amounts of neutrons and 
protons. Therefore, also the S~ hyperons appears at lower densities than the S°, here 
slightly above 3rio already. For comparison, the composition for the standard relativistic 
mean-field (RMF) model is plotted also in Fig. 3. Here, equal hyperon potentials have 
been adopted for all hyperons. The pattern of the onset of hyperon populations is 
quite similar in the two models up to 4no. For larger densities, the other hyperons are 
present at lower densities in the RMF model compared to the chiral model. However, 
the maximum density reached in neutron star configurations is about 6n , so that the 
difference between those two models will be only important for the most massive neutron 
star configurations close to the maximum mass. 

Hyperons have a dramatic effect on the nuclear equation of state (EoS), the relation 
between the pressure and the energy density of matter, which serves as the crucial input 
to the structure equations for compact stars, the Tolman-Oppcnhcimcr-Volkoff equations. 
Fig. 4 shows the equation of state for the chiral model (thick lines) and the RMF model 
(thin lines). The upper curves are calculated with nucleons and leptons only, the lower 
ones include effects from the presence of hyperons in dense matter. One sees, that 
the hyperons substantially lower the pressure for a given energy density. This effect 
is stronger than the difference between the two models used here. In particular, the 
high-density EoS with hyperons is nearly the same in both models. 

The lowering of the pressure has a destablizing effect for neutron stars, as pressure is 
needed to counteract the attractive pull of gravity. Hence, one expects that the maximum 
possible mass for compact stars with hyperons is significantly lower than the one for 
neutron stars with just nucleons and leptons. Fig. 5 demonstrates this effect of hyperons 
on the mass-radius diagram of compact stars. Matter with nucleons and leptons only 
arrive at maximum masses of 1.84M Q for the chiral and 2.16M Q for the RMF model, 
respectively. The maximum masses with hyperons included are reduced to 1.52M Q for 
the chiral and 1.55M Q for the RMF model. Note, that the radii are different already 
at lower masses which is due to the difference in the EoS slightly above normal nuclear 
saturation density. The drastic reduction of the maximum mass due to hyperons is known 
to be a quite generic feature, see e.g. [11]. 

Finally, the transition to hyperon matter is studied by tuning up the hyperon-hyperon 
interaction, which is scarcely known from the few double A hypernuclear events available. 
With increasing hyperon-hyperon attraction, the transition to hyperons becomes a first 
order phase transition. The mixed phase, with slowly rising pressure as a function of 
energy density, will make the compact star less stable. The onset of the pure hyperonic 
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Fig. 4. - The equation of state of neutron star matter in /3-equilibrium with and with out 
hyperons (taken from [8]). Thick lines show the results for the chiral model, thin lines the ones 
for the RMF model using set TM1. Hyperons substantially reduce the pressure for a given 
energy density, i.e. they soften the equation of state. 



phase with a steeply rising pressure can lead to another stable sequence of compact stars. 
For a sufficiently huge attraction between hyperons, hyperonic matter becomes more 
stable than pure nucleonic neutron star matter and the mass-radius relation changes to 
the one for selfbound stars, i.e. the mass increases from the origin with the radius as 
R 3 , so that the average density in the star is nearly constant. Fig. 6 depicts the mass- 
radius diagram for an increased hyperon-hyperon attraction. For moderate attraction, 
the mass-radius relation has two distinct branches with correspondingly two maximum 
masses. The new branch at smaller radii is another stable solution to the TOV equations, 
constituting the so called third family of compact stars [13, 14]. For even larger hyperon 
attraction, the mass-radius relation changes to the one for selfbound stars: the mass- 
radius relation starts at the origin. Selfbound stars can still possess an outer nuclear 
crust, which is determined by the low-density nuclear EoS below neutron-drip density. In 
that case, the mass-radius relation for the low-mass configurations interpolates between 
the selfbound star configuration and the normal low-mass neutron star sequence. 
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1.0 

M/M 

Fig. 5. - The mass-radius relation for the relativistic mean-field model using the parameter set 
TM1 (thin lines) and for the chiral SU(3) model (thick lines) with and without hyperons (taken 
from [8]). Hyperons substantially reduce the maximum possible mass. 



3. Strange Bosons in Hadron Stars 

Besides hyperons, other hadronic particles can be present at high density in the core 
of neutron stars. In the following, we discuss how bosons and the phenomenon of Bose- 
Einstein condensation can be described in a simple relativistic field-theoretical approach, 
here for the case of kaon condensation in neutron star matter [15, 16]. We follow the 
field-theoretical model of Ref. [17, 18]. 

In dense matter, kaons, and in particular the negatively charged K~ , can appear in 
neutron star matter. For that to happen, it must be energetically favoured to replace 
electrons by K~ meson which translates to the condition that the effective energy of 
the K~ in the dense medium must be equal to the chemical potential of electrons, 
E* K {n) — The in- medium shift of the mass and the energy of the kaons can be 
modelled by Yukawa coupling terms to scalar and vector fields generated by the nuclear 
matter. Here, for simplicity, we adopt the standard relativistic mean-field model for the 
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Fig. 6. - The mass-radius relation for different strength of the attractive hyperon-hyperon inter- 
actions as given by the coupling constant g a * (taken from [12]). A new stable solution appears 
in the mass-radius relation with similar masses but smaller radii. Selfbound hyperon stars are 
shown by long-dashed lines. 



nuclear part and just add those coupling terms for the kaon field to the Lagrangian: 
£ = '^2^>b (i^ffid^ - m B + g aB <? - g^Bl^V^ - 5 p btb-R m ) *s + -d^ad^a - ^x^o 2 

B 

(15) -U(a) - \v^V» v + \mlV^ + U(V) - + \m 2 p R^ 

where the selfinteraction terms between the meson fields read 



(16) 



1 



1 



U(a) = ^bm{g a af + -c(g„a) 4 , U(V) = -{V^f . 



If the critical condition for the onset of kaon condensation is fulfilled: 



(17) 



U K — fJ-K 
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then processes like 

(18) e~^K~ + v e n^p + R- 

produce kaons in the dense medium. The Lagrangian for the kaon field can be cast in 
the form 

(19) C K =V* il K*V> x K -m* K 2 K*K 
where the vector fields are coupled minimally 

(20) = -+- ig^KVp. + i9 P Kf K Rp 

and the effective mass of the kaon is defined as a linear shift of the mass term by the 
scalar field 

(21) m*K = m K -g aK <j . 

The in- medium effective energy of the kaon can be read off from the energy-momentum 
relation by using a plane wave ansatz for the kaon field: 

(22) lo k = m K - g aK a - g^nVo - g P KRo,o ■ 

The total energy density and pressure can be calculated with standard techniques from 
the energy-momentum tensor (assuming an ideal fluid), so that 

(23) e = ejv + ejf + e e , M 

(24) p = p N + Pe,n ■ 

The explicit expressions for the energy density originating from kaon condensation reads 
simply 

(25) e K = m* K n K . 

where uk stands for the number density of kaons. Note, that there is no direct con- 
tribution from the kaon condensate to the pressure, as it vanishes for a Bose-Einstein 
condensate. Finally, one has to fix the coupling constants for the kaons. We adopt again 
symmetry arguments to relate the vector coupling constants for the kaons to the ones of 
the nucleon 

(26) g^K = ^ff^Ar and g pK = g pN . 

Here, the coupling to the vector mesons are assumed to be universal, which is the case 
when gauging the Lagrangian, so that the coupling constants are just related by isospin 
symmetry. An extension to SU(3) is straightforward, but not necessary here, as we 
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Fig. 7. - The kaon and antikaon energy in dense nuclear matter (taken from [24]). The kaon 
energy is shifted up, while the antikaon energy is greatly reduced as a function of density due to 
an overall strongly attractive antikaon-nucleon potential. RMF: relativistic mean-field model, 
CHPT: chiral perturbation theory with different values of the £ term and for a vanishing range 
term, D — (see [24] for details) . 



discuss first neutron star matter with nucleons, leptons and kaons only and ignore the 
effects from hyperons. The scalar coupling constant can be related to the relativistic 
kaon in-medium potential via 

(27) U K (po) = -.9<xW/°o) - 9lukV {po) 

which is a more appropriate (and pragmatic) control parameter. The kaon potential 
is not well known, but likely to be considerably attractive at high densities. Coupled 
channel calculations for kaons in dense matter arrive at values of up to —120 MeV at 
normal nuclear matter density [19, 20, 21] although much shallower potentials are found 
in selfconsistent treatments, see e.g. [22, 23]. 

Fig. 7 shows the in-medium energy of kaons and antikaons as a function of baryon 
density for different approaches. The model outlined here is denoted by the label RMF. 
The prediction from the RMF approach are in good agreement with the one from chiral 
perturbation theory (ChPT), where several cases are shown but not discussed here in 
more details. The kaon energy is shifted up in the nuclear medium due to the repulsive 
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Fig. 8. - The equation of state for nucleon star matter with kaon condensation (from [18]). 
For large attractive kaon potentials, the transition to kaon condensation is of first order. The 
pressure decreases as function of energy density at the onset of kaon condensation (red dotted 
lines) indicating an instability. A Gibbs construction has to be used for the description of the 
mixed phase (blue dashed lines). 



vector potential. The antikaon on the other hand experiences an attractive vector po- 
tential, so that the effective energy decreases dramatically as a function of density. For 
comparison, a coupled channel calculation for the K~ is also plotted, with rather similar 
results to the ones of the RMF model. 

The appearance of kaon condensation in neutron star matter is accompanied by a 
strong first order phase transition. The Gibbs criteria for handling a phase transition 
with two conserved charges, baryon number and charge, states that the pressure in the 
two phases has to be equal for the same chemical potentials 



The equation of state for neutron star matter is depicted in Fig. 8. The solid line stands 
for the purely nucleonic EoS, the dotted lines for the pure kaon condensed phase with 
nucleons and kaons for different values of the kaon potential. The pure kaon condensed 
phase is seen to be unstable, there is a region where the pressure decreases with increasing 
energy density. The Gibbs construction for the mixed phase of pure nucleon matter and 
pure kaon condensed matter (dashed lines) interpolates between the two phases in a 
thermodynamical consistent way, so that the pressure is continuously rising as a function 



(28) 




,ii 



Ms = Mb » Me = Me 
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Density p (fm 3 ) 

Fig. 9. - The population of nuclear star matter with kaon condensation for a kaon optical of 
U — — 120 MeV at n — no (taken from [18]). Kaons appear around 3no and replace the electrons. 
The proton fraction balances the negative charge of the kaons and can be even larger than the 
neutron fraction for high densities. 



of energy density. 

The composition of neutron star matter with kaon condensation is shown in Fig. 9 
for a kaon potential of Uk = — 120 MeV at no. The K~ set in at 3no and reach a very 
high fraction already for slightly larger densities. As soon as a the K~ are present, the 
electron and muon fraction decreases accordingly. At 5n there arc equal fractions of 
neutrons, protons and K~ , so neutron star matter would be more aptly called nuclcon 
star matter. At even larger densities, the K~ dominate the population together with the 
protons, which ensures overall charge neutrality. The K~ in combination with protons 
are favoured in comparison to the neutrons as they are deeply bound in dense matter. 

Fig. 10 depicts the corresponding mass-radius relation by using the equations of state 
as seen in Fig. 8. The lower the kaon optical potential, the lower is the critical density 
for the onset of kaon condensation. Hence, the mass-radius relation will deviate from the 
one of the canonical neutron star at lower masses for larger values of the kaon potential. 
As kaons are in a Bose-Einstein condensate and do not give a direct contribution to the 
pressure, less mass can be supported against the gravitational pull. For low values of the 
kaon potential, the mass-radius curve simply gets unstable as soon as kaons are part of 
the composition in the core of the compact star. For the case Uk = —130 MeV, the mass 
slowly changes as a function of radius and very small radii well below 10 km are reached. 
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Fig. 10. - The mass-radius relation for kaon-condensed nucleon stars (taken from [18]). The 
maximum mass is reduced and extremely small radii of about 8 km are possible when kaon 
condensation is present. 



High compression is needed to stabilise a kaon condensate compact star. For even more 
attraction, the case Uk = —140 McV, there appears a significant difference in the mass- 
radius curves between the (thermodynamically inconsistent) Maxwell construction and 
the Gibbs construction. A strong first order phase transition is present in the core of 
the compact star from neutron star matter to kaon condensed star matter. Interestingly, 
the deviations in the mass-radius curve due to the different descriptions of the mixed 
phase remains to be just in some intermediate region of the mass-radius curve. For the 
last sequence of stable compact star configuration up to the maximum mass, the two 
descriptions give about similar results for the mass-radius curve. The unstable region in 
the case of the Maxwell construction, indicated by a change of the slope of the mass- 
radius diagram, is absent in the case of the Gibbs construction. Note, that for the Gibbs 
construction, the kaon condensed phase appears at much lower values of the mass of the 
compact star compared to the Maxwell construction. We will discuss phase transitions 
in more detail with regard to the onset of the quark matter phase at high densities, to 
which we turn now. 
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4. Strange Quark Matter in Compact Stars 

In this section, we discuss the properties of compact stars with quark matter. First, 
pure quark stars are addressed and the mass-radius diagram for so-called selfbound stars. 
Then, the quark matter equation of state is matched to the one of the low-density 
hadronic equation of state. Compact stars with both types of matter, hadronic and 
quark matter, are dubbed hybrid stars. 

4'1. General properties of quark stars. - QCD predicts that quarks at large energy 
scales are asymptotically free. For high energy-density matter one considers a gas of free 
quarks with corrections from one-gluon exchange. The pressure at zero temperature and 
finite quark-chemical potential \x reads for massless quarks 



(29) p( M ) = 



4tt 2 



(t) 



G + N f In — 

TT 



11 



2 

3*' 



In*' 
M. 



where Nf stands for the number of flavours (here Nf — 3). The constant G is scheme 
dependent and in the MS scheme given by G = G - 0.5367V/ + N f lnN f , G = 10.374 ± 
0.13 [25]. The renormalisation scale A will be fixed to be proportional to the quark- 
chemical potential. The expression for the pressure of massive quarks can be found in 
[26]. 

In compact star matter, weak reactions are in equilibrium so that 



(30) 



d 
s 

s + u 



u + e~ + u e -, 
u + e~ + v e - , 
d + u, 



and the chemical potentials for down and strange quarks must be equal, while the one 
for up quarks is just shifted by the electrochemical potential: 



(31) 



Ms = [i-d = Mm + A*e 



Note, that this amounts to baryon and charge number conservation, the chemical po- 
tentials of the quarks can be also expressed in terms of the baryochemical potential and 
the electrochemical, as done in the previous sections. This fact will be beneficial for our 
discussion of the matching of the two phases in the next subsection, for example the 
connection [i n = n u + 2^ can then be easily made. From the pressure one can compute 
the number density n and energy density e via the standard thermodynamic relations for 
each quark separately 



(32) 



d Pf(»f) 
dfXf 



and ef(nf) = -pf(fif) + fJ, f Pf{fJ, f ). 
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In addition, pure quark star matter should be charge neutral so that for the total charge 
density 

2 11 

(33) n c = -n u - -n d - -n s - n e = 

where rn stands for th number density of the species i. Finally, the total energy density 
and pressure of the quark matter is 

(34) e Q = e u + e d + e s + e e , 

(35) PQ = Pu+Pd+Ps+Pe ■ 

For comparison, one arrives at the MIT bag model by setting the coupling constant to 
zero and shifting the energy density and pressure of the quark matter by the bag constant 
B which describes the non-perturbative aspects of the QCD vacuum: 

(36) 6bag = Efree + B pbag = Pfree - B 

which can be combined to 

14 

(37) Pbag = T^bag - g# 
Phcc = e frcc/3- 

In the following we consider quark matter of massless up, down, and strange quarks. 
The current strange quark mass of about 100 MeV is smaller than the scale of the quark- 
chemical potential of 300 MeV and more. Corrections are of the order of (m s /ii) 2 and can 
be safely ignored for our purposes. For the case of three-flavour massless quarks quark 
matter consists of equal number densities of up, down and strange quarks and is charge 
neutral by itself. Hence, there are no electrons in such idealised quark matter (in fact 
the corrections from the finite strange quark mass are small). The chemical potentials 
of all three quark species are the same and just given by the quark-chemical potential. 
The quark number density and the energy density is then determined by 

/r.o\ dp 

(38) n = — — e = —p+u-n 

which fixes the quark matter equation of state in parametric form. The pressure actually 
vanishes for some critical chemical potential. Hence, there is a energy density jump from 
the vacuum to the energy density of quark matter. Pure quark matter is then stable 
at this characteristic energy density. This feature of the quark matter equation of state 
can be easily seen for the MIT bag model but is also present in the interacting model. 
The characteristic energy density is then fixed by the bag constant, as eb ag = 4_B for 
vanishing pressure. 

Fig. 11 depicts the mass-radius relation for compact stars consisting of pure quark 
matter. The solid lines show the interacting case for the perturbative QCD equation of 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




Radius (km) 

Fig. 11. - The mass-radius relation of pure quark stars in two different approaches: for the MIT 
bag model (green dashed curves) and within perturbative QCD calculation (blue solid lines). 
The general form of the mass-radius relation is remarkably similar. 



state, the dashed line the one for the MIT bag model for so-called strange stars [27, 28]. 
As the pressure vanishes at some finite value of the energy density, the quark stars are 
stabilised and have about similar densities. Hence, the mass increases with the radius of 
the quark star as R 3 starting at the origin. There is a maximum mass, as at some critical 
number of quarks the gravitational pull can not be counteracted by the Fermi pressure. 
The curves can be scaled into each other. For the MIT bag model, maximum mass and 
the corresponding radius scale as £? -1 / 2 , smaller values of the MIT bag constant give 
larger maximum masses and larger radii. The result for the perturbative QCD model are 
strikingly close to the simple MIT bag model. Here larger choices for the renormalisation 
scale result in larger maximum masses and radii. Note, that pure quark stars can only 
exist if strange quark matter is more stable than ordinary nuclear matter, so that any 
hadronic mantle is transformed to strange quark matter, which is the so-called Bodmer- 
Witten hypothesis [29, 30]. 

4'2. Hybrid stars: compact stars with quark and hadron matter. - It is more likely 
that strange quark matter is not absolutely stable, so that the quark core is surrounded 



Strangeness in Relativistic Astrophysics 



21 



by hadronic matter. The compact star is then dubbed a hybrid star. The transition 
from one phase to the other in compact star matter has to be modelled via the Gibbs 
condition for phase transitions [31] which states that in phase equilibrium the pressure 
of the two phases is the same for the same chemical potentials of the two phases 

(39) p H {fJ,B,Ve) =P0(A»B,Me)- 

Here, the index Q denotes the quantity in the quark phase, H the one in the hadronic 
phase. There are two conserved charges for compact star matter, the baryon number 
and charge, so that there are two independent chemical potentials. A Maxwell construc- 
tion could only ensure that one chemical potential is continuous throughout the phase 
transition while the other would jump discontinuously. The volume fraction of the quark 
phase 

(40) X - VQ 



Vq + Vh 

can be used to calculate the total energy density in the mixed phase 

(41) e mix =xe Q + (l-x)e/f- 

The charge neutrality is now guaranteed globally, i.e. the Gibbs construction allows for 
highly charged phases whose charges balance each other. 

The Gibbs construction can be visualised by the intersection of the two planes in 
a diagram where the pressure is plotted as a function of the two chemical potentials, 
see Fig. 12. In the plane of the pressures there is one line which indicates locally charge 
neutral matter of a single component, quark or hadron matter. The line of the intersection 
of the two pressure planes defines the mixed phase, the pressure is equal in the two phases. 
One realizes that the line of intersection is off the line of charge neutral matter for the 
quark or the hadron phase alone. The volume fraction is now chosen in such a manner 
that the total mixed phase has a vanishing charge density p and is charge neutral globally: 

(42) p m ix = X PQ + (1 - X) PH = 0. 

At the beginning of the mixed phase, highly negative charged quark matter bubbles 
appear while hadron matter is slightly positively charged. At the end of the mixed 
phase, quark matter is slightly negative charged while there are small bubbles of highly 
positive charged hadron matter. This mismatch of charges is energetically favoured, 
hadron matter is more stable for about equal numbers of protons and neutrons due to 
the asymmetry energy, so that a positive charge is advantageous. 

The corresponding composition for hybrid stars is depicted in Fig. 13 as a function 
of the radius for different choices of the MIT bag constant. For low values of the bag 
constant, the core consists of pure quark matter which can be the dominant part in the 
overall composition of the hybrid star. A sizable mixed phase exists, in particular for 
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Fig. 12. - The Gibbs phase construction for two chemical potentials by looking at the pressure 
in 3D. The surface of the hadronic pressure cuts the one of the quark matter pressure along the 
phase coexistence line where both pressures are equal. The overall charge neutrality condition 
fixes a line in the hadronic pressure area and also in the quark one. Charge neutrality is ensured 
along the mixed phase line by adjusting the volume fractions of hadron and quark matter (taken 
from [14]). 



large values of the bag constant, where the pure quark core is absent in the hybrid star. 
For sufficiently large values of the bag constant, the mixed phase as well as the pure 
quark phase is not present in the compact star, which is then an ordinary neutron star 
without any quark matter component in its core. The right plot shows the case, when the 
quarks have a quasiparticle effective mass due to interactions via gluon exchange which 
is parameterised by the coupling constant g. The onset of the mixed phase and the pure 
quark phase is shifted to lower values of the bag constant. Note, that the composition of 
a compact star with quark matter is highly sensitive to the choice of the bag constant. 
Anything between a pure quark star and an ordinary compact star of hadronic matter 
only is possible due to our present limited knowledge about the properties of quark matter 
at extreme densities. 

Finally, the mass-radius diagram for hybrid stars is shown in Fig. 14 for different 
values of the bag constant. For large radii, a mixed phase is formed in the core of the 
compact star. The mass is decreasing for lower values of the bag constant (note the small 
differences in the values for the bag constant). For some values of the bag constant, the 
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Fig. 13. - The composition of hybrid stars for different values of the MIT bag constant B (taken 
from [14]). For large values of B, the compact star is purely hadronic. A pure quark core is 
present for small values of B. Note, that for even lower values of B, the pure quark phase will 
extend up to the surface of the compact star and a pure quark star (a so called strange star) is 
formed. 



mass decreases as a function of the central energy density (for smaller radii) after reaching 
a maximum mass. This region is unstable with respect to radial oscillations. However, at 
even larger central energy densities (smaller radii) , the mass starts to increase again and 
a new sequence of stable solutions appear. The compact stars corresponding to these 
new stable sequence constitute a third family of compact stars, besides white dwarfs 
and ordinary neutron stars. They are stabilised by the presence of a pure quark phase 
which has a sufficiently stiff equation of state so as to support very dense compact star 
configurations. It is possible to have two compact stars with identical masses but with 
different radii, so-called compact star twins [13, 14]. 

In recent years, the properties of QCD at high densities have been explored in 
connection with the phenomenon of colour-superconductivity [32, 33], for reviews see 
e.g. [34, 35]. A rich phase structure exists in the QCD phase diagram at high densities 
and low temperatures right in the regime of compact star physics. In particular, several 
phase transitions might be present in compact star matter, so that the quark matter 
core of hybrid stars might contain more than one phase transition. This research field 
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R [km] 

Fig. 14. - The mass-radius relation ol hybrid stars, compact stars with hadronic and quark 
matter (taken from [14]). Configurations exists where the core consists just of the mixed phase 
(mixed core: MC) or where the core contains pure quark matter (quark core: QC). In the latter 
case, the compact star with a quark core can constitute a third family of compact stars with 
similar masses but smaller radii than the mixed core counterparts. 



is rapidly evolving and is poised for new discoveries for the properties of compact stars 
with quark matter. We refer the interested reader to the above mentioned review articles 
for further reading. 

5. Summary 

Compact stars can have a rich structure with new and exotic phases being present 
in their cores. Hyperons, heavy baryons with strangeness, can appear in neutron star 
matter. As a new fermionic degree of freedom, hyperons lower the overall pressure so that 
the maximum mass of a compact star with hyperons is substantially smaller than the 
one of a neutron star consisting of nucleons and leptons only. Also, negatively charged 
Bose-Einstein condensate could be formed via kaon condensation at high densities. As 
the Bose-Einstein condensate does not give a direct contribution to the pressure of the 
system, the stable sequence of compact stars in the mass-radius diagram simply stops 
as soon as the kaon condensed phase appears or extremely dense compact stars are 
generated with unusually small radii. Finally, at high densities quark matter can be 
present. Depending on the onset of the phase transition and the strength of the first- 



Strangeness in Relativistic Astrophysics 



25 



order phase transition from the hadronic to the quark matter phase, a third family of 
compact stars is found in the mass-radius diagram. The compact stars of this third 
family have a core of pure quark matter and smaller radii than ordinary hybrid stars. 
The physics of dense quark matter in compact star opens exciting new perspectives in 
exploring the phase diagram of QCD at high densities. The new research facility FAIR 
at GSI Darmstadt will explore this fascinating regime of strong interaction physics in 
the near future with relativistic heavy-ion collisions with bombarding energies tuned so 
as to create compressed baryonic matter at highest baryon densities. 
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